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New quark singlets T can be produced in pairs at LHC through standard QCD interactions, 
with a large cross section for masses of several hundreds of GeV. Their charged current de- 
cays TT W^bW^b, with semileptonic decay of the W pair, give a final state £vbbjj, 
as in top pair production. The mixed decay modes TT W^bHt,HtW^b W^bW^bH, 
Tf W+bZt,ZtW'b -> W+bW'bZ, with H,Z jj, yield two more jets, ivbbjjjj, or the 
same state in the latter case for Z vv. We extend previous work examining in detail the extra 
H and Z contributions and stressing one of the salient features of the TT W^bW^b signal: the 
presence of a very energetic charged lepton. We finally compare with the discovery potential in 
ivbbbb jj final states, with four b tags. 
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1. Introduction 

Fermion singlets under S\]{2)l are not present in the minimal Standard Model (SM); however, 
their existence is experimentally allowed, provided their mixing with the standard fermions is of 
order 0(0. 1) or smaller. By far, the best known example in which singlet fermions are postulated is 
the seesaw mechanism, where neutrino singlets are introduced to explain the smallness of neutrino 
masses. Quark singlets appear in some SM extensions like those with extra dimensions [1], in 
grand unified theories [2, 3] and Little Higgs models [4]. In analogy with the neutrino case, the 
existence of one 2 = 2/3 singlet T (or more) could provide a clue to understand the largeness of 
the top mass [5,6]. 

The possible existence of a 2 = 2/3 singlet T with a mass of several hundreds of GeV would 
have very interesting phenomenological consequences. Its mixing with the top and up or charm 
quarks would generate flavour-changing neutral coupUngs at the tree level, large enough to be 
detectable at the Large Hadron Collider (LHC) [7,8]. In low energy processes, its contribution 
could lead to significant departures from SM predictions for several observables, like rare kaon 
decays, Z)" — D^ mixing and CP asymmetries in B oscillations [8-11], compatible with present 
experimental measurements. But the most interesting "indirect effect" of a quark singlet is Higgs 
boson discovery [12, 13]: in the presence of a heavy singlet, the luminosity required for Higgs 
discovery at LHC is dramatically reduced. For instance, for tut = 500 GeV, the luminosity needed 
for discovery of a 115 GeV Higgs is only 6 fb~^ [14], ten times smaller than within the minimal 
SM. 

Besides their indirect effects, Q = 2/3 singlets can be directly observed at LHC. They can be 
produced in pairs by QCD interactions [3, 13], or in association with a light jet, through its mixing 
with the bottom quark Vjb [15, 16]. Both processes have their own advantages and weaknesses. 
Pair production has a model-independent cross section determined by QCD interactions, which 
however decreases quickly with mj, due to phase space suppression. Single T production cross 
section decreases more slowly with the heavy quark mass, but on the other hand it is proportional 
to |Vri)P, and gets suppressed for small mixings. Branching ratios for Tf final states involve two 
heavy quark decays and are thus smaller than for T j production. But they also have more fermions 
in the final state, and then their backgrounds are less important. A heavy quark T can decay via 
charged currents T ^W^b, neutral currents T ^ Zt and, if the Higgs boson is light enough, 
T ^ Ht as well. The decay widths are approximately Br(r W^b) = 0.5, Br(r Zt) = 0.25, 
Br(r Ht) = 0.25 (these figures are exact for » Mw,Mz,Mh)- We neglect decays to first 
and second generation quarks, involving mixing angles Vts, etc. which are bound to be much 
smaller than Vrb [8]. The leading decay T W~^b i^vb produces a very energetic charged 
lepton, which can be used to separate the T signals from their backgrounds. The decay mode 
T ^Zt^i+i-i+vb, £J' = e ,11, gives a cleaner final state than T ^W^b, but with a branching 
ratio 10 times smaller. The remaining channel T ^ Ht ^ bbi^vb yields 3 b quarks, which can be 
tagged to reduce backgrounds. Additionally, it constitutes an important source of Higgs bosons for 
moderate mt values, for which the Tf production cross section is large. 

Single T production at LHC has been studied in detail, in the three decay modes of the heavy 
quark [17, 18]. The best results have been found for the T W~^b chaimel. For pair production, 
the modes TT W^bW-b and TT W^bHt,HtW-b have been analysed [14, 19]. Decays to 
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W^bW~b, with subsequent semileptonic decay of the W pair, give an experimental signature of 
one charged lepton, two b jets and at least two additional jets. Decays to W~^bHt,HtW~b include 
two more jets in the final state from the Higgs decay, which mostly correspond to hh pairs. The 
same occurs in TT W^bZt,ZtW^b with hadronic Z decay, but with a much smaller fraction 
of bb pairs. In the latter mode, the Z boson can also decay invisibly, giving a further ivbbjj 
contribution. In the following we will carefully analyse the additional H and Z signals, extending 
what was presented in Ref. [19], and show how they enhance the sensitivity in the W~^bW~b 
channel. 

2. Signal and background simulation 

The main backgrounds for the W~^bW~b signal and additional contributions 

gg,qq ^ TT ^W^bW'b^tvbqq'b , i = e,^, 

gg,qq TT ^W+bHt,HtW^b^W+bW^bH ^i+vbqq'bbb/cc, 

gg,qq Tf ^W+bZf,ZtW-b ^W^bW-bZ ^ i+vbqq'b (/'q" /VV , (2.1) 

are given by tt, Wbb jj, Zbb jj and tb j production. The charge conjugate processes are understood 
to be summed in all cases. In our evaluations we do not consider final states with T leptons (which 
can decay leptonically z eViVe, T — > AiVrV^) which are estimated to increase cross sections 
3% or less [14]. Wjjjj and Zjjjj production seem to be negligible after the requirement of 
two b tags, which suppresses their cross sections by a factor ~ 10^^. The signal and the tt, tbj 
backgrounds are evaluated with our own Monte Carlo generators, including all finite width and 
spin effects. We calculate the matrix elements using HELAS [20] with running coupling constants 
evaluated at the scale of the heavy quark T or t. Wbbjj and Zbbjj are calculated with ALP GEN 
[21]. The bottom quark mass nih = 4.8 GeV is kept in all cases, and we take Mh = 115 GeV. We 
use structure functions CTEQ5L [22], with = s for TT, tt and tbj, and = + Pj^^ for 
Wbbjj, Zbbjj, being the partonic centre of mass energy, and p^^^ transverse momentum 
of the gauge boson. Events are produced without kinematical cuts at the generator level in the 
case of TT, tt, while for the other processes we set some loose cuts. For tb j we only require 
pseudorapidities | ?] | < 3 for i?, j. For Wbb jj we set Pt>l5 GeV and 1 7] | < 3 for the charged lepton, 
the b quarks and the jets, and lego-plot separations ARjj,ARhj,ARhh > 0.4, ARih,ARij > 0.2. For 
Zbb jj we require p, > 15 GeV, \ri \ < 3 for b quarks and jets, |t7 | < 10 for the charged leptons, and 
ARjj,ARbj,ARbb>OA. 

The events are passed through P YTH I A 6.228 [23] as external processes to perform hadro- 
nisation and include initial and final state radiation (ISR, FSR). We use the standard P YTH I A set- 
tings except for b fragmentation, in which we use the Peterson parameterisation with Sh = —0.0035 
[24]. A fast detector simulation AT LEAST 2.60 [25], with standard settings, is used for the 
modelling of the ATLAS detector We reconstruct jets using a standard cone algorithm with 
AR = y^ATTp + ^A^)^ = 0.4, where rj is the pseudorapidity and (j) the azimuthal angle. The 
package ATLFASTB is used to recalibrate jet energies and perform b tagging, for which we select 
efficiencies of 60%, 50% for the low and high luminosity LHC phases, respectively. Efficiencies 
for charged lepton identification and triggering are not included. The simulated events are required 
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to fulfill these two criteria: (a) the presence of one (and only one) isolated charged lepton, which 
must have transverse momentum pt > 20 GeV and < 2.5; (b) at least four jets with pt > 20 
GeV, |t]| < 2.5, with exactly two b tags. The presence of the high-/?, charged lepton provides a 
trigger for the events. The cross sections times efficiency of the five processes after these pre- 
selection cuts are collected in Table |l|, using a 60% b tagging rate. We observe that the additional 
H and Z contributions to the signal actually have a larger cross section than the signal itself. For 
example, in the first case of mr = 500 GeV we have Br(r W+b) = 0.503, Br(r Ht) = 0.331, 
Br(r Zt) = 0.166. The W^bW^b mode has a branching ratio of 0.253 (without including W 
decays), while the mixed modes have branching ratios of 0.333 and 0. 167, respectively, which have 
to be multiplied by extra factors Bi{H jj) = 0.96 and Br(Z jj, vv) = 0.90, respectively. 

Process a x eff 



rr (500) 37.3 fb + 46.5 fb(//)+ 19.8 fb(Z) 

Tf (1000) 0.618 fb + 0.638 fb (H) + 0.481 fb (Z) 

tt 18.8 pb 

Wbbjj 1.23 pb 

Zbbjj 246 fb 

tbj 710 fb 

Table 1: Cross sections of the TT signals (with nij ~ 500, 1000 GeV) and backgrounds after pre-selection 
cuts. 



The TT signal can be discovered by the presence of peaks in the invariant mass distributions 
corresponding to the two decaying quarks. In order to reconstruct their momenta we first identify 
the two jets j\, 72 from the W decaying hadronically. The first one j\ is chosen to be the highest 
Pt non-b jet, and the second one 72 as the non-^ jet having with j[ an invariant mass closest to 
Mw- The missing transverse momentum is assigned to the undetected neutrino, and its longitudi- 
nal momentum and energy are found requiring that the invariant mass of the charged lepton and 
the neutrino is the W mass, {p£ + Pv^ = M^. This equation yields two possible solutions. In 
addition, there are two different pairings of the two b jets to the W bosons decaying hadronically 
and leptonically, giving four possibilities for the reconstruction of the heavy quark momenta. We 
select the one giving closest invariant masses mj?'', m^^'' for the quarks decaying hadronically and 
semileptonically. Obviously, this reconstruction procedure does not require the knowledge of the 
heavy quark mass. 

3. Results for mr = 500 GeV 

We first present our results for a heavy quark mass of 500 GeV. The reconstructed masses of the 
T quarks decaying hadronically and semileptonically are plotted in Fig. |l| (a) and (b), respectively, 
for the W^bW^b signal, the H and Z contributions and the tt background. All histograms are 
normalised to 2000 events. The mass reconstruction works reasonably well for TT,tt W^bW^b 
decays, and can be used to separate them. For the H and Z signals the distributions are much 
broader, therefore eliminating the tt background also suppresses a large fraction of these signals. 
The same happens with the charged lepton transverse momentum p^^^, shown in Fig. [l] (c). The 
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distribution for the W^bW^b signal has a long tail, not only due to the higher energy available 
but also because the fraction of W bosons with negative helicity produced in T decays is smaller 
than for the top quark, and then the charged lepton from W decay is produced more towards the 
W flight direction (see section The momentum of the charged lepton can be used to suppress tt 
production, but unfortunately requiring high p\^^ also reduces the H and Z contributions, in which 
the charged lepton comes from t Wb decays half of the time. Our selection cuts to suppress 
backgrounds are 

^;,max ^250 GcV , > 1 50 GcV , 

pf"^ > 50 GeV , 50 GeV < p{< 600 GeV , 

//f> 1000GeV, (3.1) 

where /j/'™^" and /jf '™^" are the transverse momenta of the fastest jet and fastest jet, respectively; 
p/ is the missing transverse momentum and Hf the total transverse energy. The requirements on 
these additional variables do not significantly affect the H and Z contributions. With these selection 
cuts in mind, for our high-statistics evaluations we require at the generator level the presence of 
a charged lepton with pt > 30 GeV, a jet with pt > 200 GeV and, for Wbbjj and Zbbjj, one b 
quark with pt> 100 GeV. This last cut does not bias the sample because in these two processes the 
two non-b jets mostly originate from light quarks and gluons, for which the b mistag probability is 
very low. Thus, the Z^-tagged jets correspond to the b quarks most of the time. For the backgrounds 
we simulate the number of events corresponding to 10 fb^^ of luminosity, and for the signals 
we simulate 100 fb^^ and divide the result by 10, to reduce statistical fluctuations. Signal and 
background cross sections are scaled by K factors which approximately take into account higher 
order processes with extra jet radiation (see Ref. [19] for details). 

The mass distributions after selection cuts of the SM background and background plus signal 
are shown in Fig. The number of events N^ut corresponding to each process can be read in 
Table ^. The second column displays the numbers of events in the peak region, defined as 

300 GeV < mf"^ < 660 GeV , 300 GeV < m'^^ < 660 GeV . (3.2) 

In this region, the statistical significance of the signals is J^o = S/^/B = 10.94 or, if we include 
a 5% background normalisation uncertainty, y=, = S/y^B+ (0.05B)^ = 7.69. Discovery of the 
new quark would be possible with 2.1 fb^' of luminosity (2.8 fb^^ if we include the background 
uncertainty). 

4. Results for mj = 1 TeV 

The simulation in this case is done for 300 fb^^ in the high luminosity LHC run. The b 
tagging efficiency assumed is of 50%. The reconstructed masses of the two quarks are shown in 
Fig. ^ (a) and (b). For the H and Z signals the events are spread over a wide range of rrij^'^ and 
m'^^ values, and for this T mass a smaller fraction of events is in the vicinity of mr- The charged 
lepton transverse momentum is plotted in Fig. |3| (c). For a heavy quark mass of 1 TeV, the lepton 
momentum distribution from T W^b — > £^vb has a very long tail. Requiring p^^'^ > 200 GeV 
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(a) 



(b) 




Figure 1: Reconstructed masses of the heavy quarks decaying hadronically (a) and semileptonically (b), and 
transverse momentum of the charged lepton (c), after pre-selection cuts. 



> 



SM 

SM + r (500) 




(a) 




(b) 



Figure 2: Reconstructed masses of the heavy quarks decaying hadronically (a) and semileptonically (b), 
after selection cuts. The dashed hnes correspond to the SM predictions, while the full lines represent the SM 
plus a new 500 GeV quark. 
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Process 


Ncut 


A^peak 


TT (500) 


201.7 


135.5 


Tf (500, H) 


139.4 


59.5 


TT (500, Z) 


58.5 


26.5 


tt 


1609 


300 


Wbb jj 


287 


81 


Zbbjj 


39 


13 


tbj 


70 


16 



Table 2: Number of events A^cut passing the selection criteria and number of events A^peak passing the selec- 
tion cuts which are in the peak regions. 

eliminates 98% of the « background, while keeping 51% of the W^bW^b signal. Our selection 
criteria are in this case 

pI'p > 200 GeV , 50 GeV < i/,< 400 GeV , 

//, > 1800GeV. (4.1) 

The mass distributions after these cuts are shown in Fig. 0. The number of events A^cut correspond- 
ing to each process can be read in the first column of Table |3|, and the number of events in the peak 
region 

800 GeV < mf'^ < 1200 GeV , 800 GeV < m'^P < 1200 GeV (4.2) 

in the second column. The statistical significance of the T signal is J^q = 9.09, ^5 = 8.82. The 
luminosity required to discover a 1 TeV quark in pair production is 90 fb^ ' . 



Process 


Ncut 


A'peak 


TT (1000) 


58.2 


33.5 


TT (1000, H) 


39.6 


7.8 


TT (1000, Z) 


21.0 


5.1 


tt 


208 


10 


Wbb jj 


132 


15 


Zbbjj 


19 


1 


tbj 


3 






Table 3: Number of events A^^ut passing the selection criteria and number of events A^peak passing the selec- 
tion cuts which are in the peak regions. 



5. Discussion 

We can compare the results obtained here with those from TT —>■ W^bHt,HtW^b decays in 
the final state ivbbbbjj, requiring four b tags [14]. For nij = 500 GeV and 10 fb^' of luminosity. 
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Figure 3: Reconstructed masses of the heavy quarks decaying hadronically (a) and semileptonically (b), and 
transverse momentum of the charged lepton (c), after pre-selection cuts. 




(a) (b) 

Figure 4: Reconstructed masses of the heavy quarks decaying hadronically (a) and semileptonically (b), 
after selection cuts. The dashed hnes correspond to the SM predictions, while the full lines represent the SM 
plus a new 1 TeV quark. 
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the latter channel has a similar sensitivity, J^o = 10.86 , = 10.13, which is achieved with a 
good b tagging performance, which suppresses the dangerous ttjj and W + jets backgrounds. We 
also note that the ivbbbbjj channel has a slightly larger branching ratio than ivbbjj final states. 
Combining both channels, a 500 GeV singlet can be discovered already with a luminosity of 1.2 
fb^^ 

For niT = 1 TeV the situation is rather different. The TT W^bW^b signal is characterised 
by the presence of a very energetic charged lepton, not only due to the higher energy available in T 
decays but also because of spin effects. Let us analyse these issues more in detail. For mj- = 1 TeV, 
the fraction of W bosons with negative helicity produced in T decays is /_ = 0.013, compared to 
/_ = 0.297 for the top quark, and then the charged lepton from W decay is produced more towards 
the W flight direction. (The fraction of W bosons with positive helicity /+ is negligible in both 
cases due to angular momentum conservation; the zero helicity fraction is then /o ~ 1 — /_.) This 
is reflected in the charged lepton energy distribution in the rest frame of the decaying quark, which 
is [26] 

1 dV 1 

r dEf 



max 

i 



.^mmp 



7min\2 



U + 3iEr-Eeff- 



+6{Er''-Ei){E,-Er)fo\ 



(5.1) 



with f™™, EJ^™ the maximum and minimum charged lepton energies in the rest frame of the 
decaying quark {Ef'' = 18.5 GeV, E^^^" = 87.4 GeV for t; = 3.2 GeV, 5^"^^" = 500 GeV for 
T). The normalised energy distributions of the charged lepton from t and T semileptonic decays are 
plotted in Fig. |5[ where we have scaled both axes to fully appreciate the spin effect. The different 
shape of the curves is reflected in the p^^^ distributions in Fig. ^. 




0.2 0.4 0.6 



Figure 5: Normalised energy distributions of the charged lepton from t and T semileptonic decays. 



For mj = 1 TeV, the requirement of high p, can be used to reduce backgrounds significantly 
in the Tf W^bW^b mode but not in TT —>■ W^bHt,HtW^b decays, where half of the events 
have a charged lepton coming from a top (anti)quark. (By the same reason, the H and Z contribu- 
tions studied in this work are suppressed.) This might be partly compensated with a good b tagging 
efficiency to reject backgrounds, but at the high luminosity LHC run, where masses of the order of 
1 TeV can be investigated, light jet and charm rejection is degraded. A preliminary analysis has 



9 



New signals in pair production of heavy 2 = 2/3 singlets at LHC 



J. A. Aguilar-Saavedra 



shown that for mr = 1 TeV and 300 fb of luminosity the significance achieved in the Ivbbbb jj 
channel is S/s/B < 3, three times smaller than in the ivbb jj mode analysed here. 
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